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The list of textbook tunneling formulas is extended by deriving exact expressions for the transmission coefficient in 
graphene ribbons with armchair edges and the step-like and barrier-like profiles of site energies along the ribbon. These 
expressions are obtained by matching wave functions at the interfaces between the regions, where quasiparticles have 
constant but different potential energies. It is shown that for an Uq high barrier and low-energy electrons and holes, 
the mode transmission of charge carriers in this type of ribbons is described by the textbook formula, where the con- 
stant barrier is replaced by an effective, energy-dependent barrier, Uq — * U(E). For the lowest/highest electron/hole 
mode, U (E) goes, respectively, to zero and nonzero value in metallic and semiconducting ribbons. This and other 
peculiarities of through-barrier/step transmission in graphene are discussed and compared with related earlier results. 
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1 Introduction. In recent literature, a considerable at- 
tention has been paid to modeling of charge transport in 
graphene HE ,34,5,678,9 Jolll II . This report gets in fo- 
cus two classic problems of Quantum Mechanics which 
have been discussed in Refs. [3.4,5,6]: Particle transmis- 
sion in a potential that has a step-like or barrier-like pro- 
file. The transmission coefficient has been shown essen- 
tially different from the textbook formula I12H . In partic- 
ular, the normal incidence of electrons or holes on the in- 
terface between two regions with different potential ener- 
gies results in the full transmission without backscattering 
J3H. In the Dirac relativistic quantum mechanics, a sim- 
ilar behavior of the transmission of massless fermions is 
known as the Klein paradox |[T3l . Because the quoted ana- 
lytical treatments were based on the use of the Dirac equa- 
tion and matching either the wave functions or transfer ma- 
trices at the interfaces of different regions, the rederivation 
of the Klein result was natural rather than surprising. As 
is well known, this approach for graphene is restricted to 
energies around the point of neutrality that implies that the 
long wave limit is valid 1141 . However, it is not at all en- 
sured a priori that all calculations in this approximation 
lead to the results which follow from exact calculations and 



then, passing to the long wave limit. Here, the transmission 
coefficient is found by exploiting the matching technique, 
as in Refs. 1U[5]|6), but taking into account the exact en- 
ergy spectrum of graphene. In other words, we solve the 
Schrodinger equation to which the Dirac equation is a cer- 
tain approximation. Our general conclusion is that the va- 
lidity of relativistic approach to the description of charge 
transport in graphene ribbons is limited, strictly speaking, 
to the zero-mode transmission in metallic ribbons. 

Our model system is an armchair graphene ribbon (GR) 
with site energies taking zero and Uq values as illustrated in 
Fig. 1. This choice is dictated by the following considera- 
tions. First, armchair GRs can have either metallic or semi- 
conducting spectrum depending on the ribbon width lfT31 . 
Thus, the quantum conductance of both basic graphene 
materials can be studied on equal footings. Second, distinct 
from zigzag GRs, the armchair GR spectrum does not have 
a special band of edge states that complicates the descrip- 
tion. Third, the long wave limit of the transmission coef- 
ficient (obtained here with the account to the discreteness 
of graphene ribbons) can be compared with earlier derived 
formulas that makes this analysis particularly instructive. 
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Figure 1 Upper part: An armchair graphene ribbon with a step-like (left) and barrier-like (right) profiles of site energies; 
gray and black circles corresponds to C atoms with site energies equal to zero and Uq, respectively (filled by hydrogens 
dangling bonds along edges are not shown). Blue (red) shaded regions indicate a higher (lower) concentration of electrons. 
Lower part: Schematic representation of the singlemode and multimode electron transmission in an armchair graphene 
ribbon. A step-like and barrier-like change of site energies Uq models a gate voltage V g , Uq = eV g ; red and black horizontal 
lines correspond to the neutrality point, E = Ep = 0. On the right, only conduction bands are shown. The potential 
difference associated with a voltage source that drives electrons from the left to the right is not shown, a « 0.246 nm. 



2 Transmission coefficient. The solution to the sta- 
tionary Schrodinger equation H'P = E'F with the nearest- 
neighbor tight-binding Hamiltonian for a honeycomb lat- 
tice has the form (C-C hopping integral \t\ is a unit of en- 
ergy) 



(1) 



m— 1 n— — oo a 



where |m, n, a) is the 2p z orbital at the ath atom of ben- 
zene ring with coordinates {to, n}, a = l,r,X, p with the 
meaning of labels explained in Fig. 1; |A/"+l, n, a=l, r) = 
0. According 021, coefficients ^ i „. q=ZiI .= sin(£j-m)0j 
& = nj/{Af+l), j = 1,2,.. .,AA, and ' 
(the site energy is zero) 



n , a ' 



m,n,a— I,', 

J n a satisfies equation 



bn,a — ffa,;^n-l,r + .9a, r^n+l^p 



(2) 



where gf r 



rtj, six 



9 J r,r, Vj9i, r = 4c0S 2 (£ J /2), 

E[E 2 - 1 - 4cos 2 (^/2)], and zeros of V 3 = 
[E 2 — 4cos 2 (£j/2)] 2 — E 2 determine the n electron spec- 
trum of an A^-long acene, C^A/^tbAf-M- For site energies 
equal to Uq, E E = E - Uq. 

As shown in I16II17I . states of ir electrons in armchair 
ribbons can be classified in 2Af "j-minus" and "j-plus" 



conduction ID bands E^(n^~) and equal number of va- 
lence bands —E^(k^~). Since we are interested in the wave- 



like solutions to the above equation, 
where satisfies the dispersion relation 



cf)^ exp(iK ■ n), 



E 



±2 



1 ±4cos(^-/2)cos(k^/2) + 4cos 2 (£j/2), (3) 



it is convenient to introduce the phase shift between the 



Zth and ?'th sites: A 



-4) 



l - Yr- - e Wj = Sf,i/(1 - 
9r )■ ^ W ^ be seen soon that the change of this phase, 
6j—>9j, that corresponds to a change of site energies by 
Uq, E^^E^ — Uq = E^ and Kf—>tif, plays an impor- 
tant role in determining electron transmission to/through a 
region, where the site energy equals Uo. 

With reference to Fig. 1, the wave function fl}, describ- 
ing incoming from the left and reflected backward or trans- 
mitted to the right electrons, can be represented as 



e iK f n + rj e- iK t n , n <l, 
Ue iH t n , n>l, 



(4) 



9 , 

n,l 



e i[it^(n-l)+^] +r . e -i[K J ± (n-l)+e j ]^ n< lj 



(5) 
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for the step-like potential, and 



in^n i „ — ire- n 

e i +r 7 e ' , 



n<l, 



t Je lK f n , 



rK f n ,n€l,N, (6) 
n>iV, 



aj . e *[«*(n- 1 )+^] 

^. e i[s ± («-i)+e J ] 



6,- e ;^ ±( "- 1)+e1 ,n G T^, 



n>7V 



(7) 

for the barrier-like potential. Matching the wave functions 
at the interfaces between regions with different site ener- 
gies gives us the necessary equations for unknown coeffi- 
cients dj, bj, rj, and tj. Finding the amplitude of transmit- 
ted j wave yields the transmission coefficient for the given 
mode, 



Esin{Rf/2)/[E sin(«f /2)] | , 



1. 



(8) 



where the upper and lower lines refer to the potential step 
and barrier, respectively. The total transmission coefficient, 
T(E) = Tj is determined by the number of "open" 
modes, see below. 

2.1 Step U . By exploiting tj from Eqs. © and ©, 
we obtain 

I sin 9j sin C 



T, 



" sin 2 [(^±^)/2]' 
or, rewritten as a function of wave vectors, 



cos 2 (£j/2) sm(nf/2) sin(^/2) 
cos 2 fe-/2)sin 2 [( K ±±K±)/4] - (£/o/4) 2 



(9) 



(10) 



The upper (lower) sign in these equations corresponds to 
E > (E < 0); energy is supposed to be positive; due to 
the spectrum symmetry Tj{E, Uq) = Tj(—E, — Uq). 

Note that dispersion relation Ej (R,j ) can be satisfied 
by both real and imaginary values of kj. However, ac- 
cording to Eq. © (and in analogy with the textbook treat- 
ment 11211 ) imaginarity of kJ~ makes Tj zero. Distinct from 
the textbook case is that the unit transmission can occur 
at 6j = dj ; fulfillment of this equality does not necessar- 
ily requires = that is the usual condition of the unit 
transmission. 

Here, our prime interest concerns the energy region 
close to the Fermi energy Ep — of undoped graphene, 
where the long wave approximation, kj~, \E\ « 1 (also 
implying \Uq\ << 1), provides a reliable description. For 
this energy region, the energy scale \/3|i|/2 is more con- 
venient. Henceforth, it is used instead of \t\ together with 
a new notation V3k x = K,J—j,± IJ with [i = 0,1,... << TV. 
The new variable satisfies the following set of dispersion 
relations, 



E = ±\ k 



Af- 



_ j |±m| s i = Mi 

i = 1,2,4,5, 



, 3/x— 1, 3^+1, . 



(11) 



which are equivalent to Eq. (© in the long-wave limit lfl8l . 
The upper line of Eq. dTTb refers to metallic GRs, j* = 
2(A/"+l)/3 is an integer; lower line refers to semiconduct- 
ing GRs, where j* = (27V + 3) /3 or j* = (27V + l)/3 are 
integers. 

For small energies, Eq. ( TTOb takes the form, 



T 



u 2 \ 



I {k x ± k x ) 2 

which can be rewritten with the use of Eq. (fTTb as 



Ti 



±y/(E 2 -k 2 )(E 2 -kf) +EE- k\ 



(12) 



(13) 



According to this equation, the zero-mode transmission in 
metallic GRs has the unit probability, To = 1, independent 
of the value of Uq. This specifies the absence of backseat - 
tering under the normal incidence in graphene n-p junc- 
tions [3 1 in the context of metallic graphene ribbons. 

The obtained result is worthwhile comparing with the 
textbook formula for the probability of over-step transmis- 
sion, D = VEE/ [VE + VE^j lIHl. The difference be- 
tween transmission in metallic and semiconducting GRs 
is substantial only in the case of singlemode or few mode 
transmission. For energies E»\Uq\, that is E xiE, Ti «1 
for the mode majority, independent of whether the ribbon 
has a metallic or semiconducting spectrum. The general 
behavior of T=£)fT Ti, i max »1, as a function of U 
is as follows: For \U \ = qE, q » 1, r~2/g; T(U =0)» 
2 -max. T(U =E)=0; and T([/ =2£;)=T max <2i max . 

2.2 Barrier U . By finding tj from Eqs. ©, (01 and 
substituting it in Eq. $Q we obtain 



sin 



sm 



sm 



+ ( 



cos f 



COS 



9jY S m 2 (RfN) 



(14) 

Similarly to Eq. (0, this representation is characteristic for 
graphene structures. It shows, in particular, that the unit 
transmission occurs under the coincidence of phases 6j 
and 9j . This can be regarded as a sort of new resonances 
which differ from the familiar condition of resonances for 
the over-barrier transmission, sin(S^TV) = 0. 

In terms of wave vectors, Eq. (fl4T i has the form 



Ti = 



• 2 ± • 2 -± 

sin Kj sin Kj 



sin 2 ref sin 2 Rf + Uf sm 2 (nf N) ' 



(15) 



where (in units of \t\) 



2 



[EE- 1+4 cos 2 (&/2)] 



cos(«f /2)cos(«f/2) 



cos 2 (6/2) 



(16) 
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To facilitate the comparison with earlier results, it is 
instructive to look at these expressions for small energies 
of electrons and holes. In the long-wave limit, an analogue 
of Eq. ( fT~2l > for the barrier-like potential reads 



Ti = 



1,2 U2 



kfU§ sin 2 (V3k x N)' 



(17) 



Again, To = 1 for metallic ribbons. Except this case, under 
the replacement fcjf/o — ► Uq/2 (and the above mentioned 
convention regarding energy units) the above equation co- 
incides with the textbook formula for the through/over- 
barrier transmission [12 1 - 

If the energy of incident electrons is tuned to the neu- 
trality point of the scattering region, E — Uq, it follows 
from Eq. ( TTTb 



Ti = 



k 2 x + Ul sinb 2 (VZhN) ' 



(18) 



Formula (4) derived in Ref. J6) for this case reads (in orig- 
inal notations) 



1 



cosh 2 Lq n + (<7„/fcoo) 2 sinb 2 Lq r 



(19) 



where, to our understanding, k 2 ^ = k 2 + g 2 , with k x and 
q n = ir(n + l/2)/W, n = 0, 1, having the meaning 
of the longitudinal (along the ribbon) and transverse com- 
ponents of the wave vector, respectively. Whatever reading 
used, we could not agree our Eqs. ( fT7] i and Eq. (TT~8b neither 
with the above equation nor with Eq. (8) from the same 
reference. 

As seen from Eqs. (T% and ( fT8l . the multimode and 
singlemode transmissions must be distinguished. In the mul- 
timode transmission, E ps fc,m»x, i max >> 1, and for 
\E\ » \U \, Ti(E) ps 1 for all open modes, \E\ > k t . 
In this case, the total transmission coefficient for metallic 
and semiconducting GRs differs only marginally. In con- 
trast, in the case of singlemode transmission, the differ- 
ence is substantial. Because of the zero and finite values of 
ko, T = 1 for metallic GRs, but for semiconducting GRs, 
T ~ exp(-^), if fc < \E\ « fei and L/W > 1. 
The latter result agrees with the expression for the Green 
function of N x Af honeycomb lattice lfl9l . 

Another formula worth mentioning in the present con- 
text refers to an expression for the tunneling probability 
that was obtained in [41 for massless Dirac fermions with 



(dimensional) kinetic energy E = ±.TiVF^jk 2 + k 2 in a 

2D space — oo < x, y < oo, where the potential energy 
is Uo, if < x < L, and zero otherwise. Assuming the 
equivalence of this expression for energy with Eq. (fTTT i and 
by using the correspondence \/3a\t\/2 — Hvf and 



it can be proved that Eq. ( TPTI i and an expression for T 
that follows from Eq. (3) in Ref. J4| have exactly the same 
form. 

To wind up this report, we would like to emphasize that 
the methodology used for the derivation of the new formu- 
las for the transmission coefficient in armchair graphene 
ribbons is the same as in J3]|4]|5]|6l ■ Distinct from the pre- 
vious considerations is the use of exact solutions of the 
model Hamiltonian for the description of tt electron states 
in ideal armchair ribbons. 
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